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Abstract We address the problem of the Euclidean upgrad- 
ing of a projective calibration of a minimal set of cameras 
with known pixel shape and otherwise arbitrarily varying 
intrinsic and extrinsic parameters. To this purpose, we in- 
troduce as our basic geometric tool the six-line conic variety 
(SLCV), consisting in the set of planes intersecting six given 
lines of 3D space in points of a conic. We show that the set of 
solutions of the Euclidean upgrading problem for three cam- 
eras with known pixel shape can be parameterized in a com- 
putationally efficient way. As a consequence, we propose an 
algorithm that performs a Euclidean upgrading with 5 (the- 
oretical minimum) or more cameras with the knowledge of 
the pixel shape as the only constraint. We provide experi- 
ments with real images showing the good performance of 
the technique. 

Keywords Camera autocalibration • Varying parameters • 
Square pixels • Three-dimensional reconstruction • Absolute 
Conic • Six Line Conic Variety 



1 Introduction 

Three-dimensional reconstructions from images are often 
obtained with calibrated cameras, i.e., cameras whose pa- 
rameters have been previously computed using calibration 
objects in a controlled environment [10, p. 201]. Unfortu- 
nately, in many cases such conditions are not available, e.g., 
when the images have been acquired with non-specialized 
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equipment or taken with a different initial purpose. To obtain 
3D reconstructions without a priori knowledge of the cam- 
era parameters or the scene content, we need autocalibration 
algorithms, which allow for free variation of the camera pa- 
rameters and provide a significant amount of flexibility in a 
wide variety of scenarios: we recall that just zooming a cam- 
era varies not only the focal length but also the location of 
the principal point |fT9l . 

The seminal paper fTTl proposed the first autocalibra- 
tion method, requiring images to be taken with cameras with 
constant intrinsic parameters. A common framework for au- 
tocalibration was later provided by the concept of geomet- 
ric stratification |^^. This technique splits the camera cal- 
ibration and scene reconstruction into three steps: first, the 
recovering of a projective reconstruction, i.e., a 3D scene 
and a set of cameras differing from the actual ones in a 
spatial homography. Second, an affine reconstruction, i.e., 
a reconstruction differing from the actual 3D scene in an 
affine transformation, is obtained by finding the location of 
the plane at infinity |12|. Finally, a Euclidean reconstruc- 
tion, differing from the actual scene in a similarity, is found 
through the location of the absolute conic at infinity or any 
equivalent geometric object. A general reference for the sub- 
ject is 161, where an extensive bibliography can be found. 

In order to upgrade a projective reconstruction to an affine 
or a Euclidean one in the absence of knowledge about the 
scene, some data about the camera parameters must be avail- 
able. For example, in the originally addressed autocalibra- 
tion problem, this additional piece of data was the constancy 
of the camera intrinsic parameters [11], resulting, for each 
camera pair, in a couple of second-degree equations in the 
camera parameters known as Kruppa equations. The insta- 
bility problems of these equations have been studied in ifTSl . 

Another constraint that has been studied is that in which 
the principal point is assumed to be known. Then the dual 
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AQC 1 SLCV 


Common features 




Basic geometric object 


Isotropic lines through the optical centers of the cameras 


Intrinsic parameter assumption 


Known pixel shape (skew and aspect ratio) 


Differences 




Type of algorithm 


Linear 
(solution of a homogeneous system) 


Non-linear 
(bidimensional search using second degree equations) 


Required number of cameras 


> 10 


>5 


Optimal w.r.t. number of cameras 


No 


Yes 


Geometric object used 


Quadric of P^ 


Algebraic surface of P-** 


Degree of the geometric object 


2 


5 


Geometrical meaning 


Lines intersecting the absolute conic 


Planes with conies intersecting all the isotropic lines 


Integration with scene knowledge 


Pairs of orthogonal lines 


Parallel lines. Points at infinity 



Table 1 Comparison of the Absolute Quadratic Complex (AQC) and the Six-Line Conic Variety (SLCV) approaches. 



absolute quadric f2T\ can be found by solving a set of ho- 
mogeneous linear equations 1161 . 

In this paper we address the problem of autocalibration 
in its less restrictive setting in practice: cameras with arbi- 
trarily varying parameters with the exception of the pixel 
shape, which is assumed to be known. It can be easily seen 
(see section 2) that this is equivalent to having cameras with 
square pixels. The possibility of employing this constraint 
was studied in [7J, where an algorithm based on bundle ad- 
justment was considered. Algorithms based on this restric- 
tion have also been proposed [13, 22,231 that result in a set 
of linear equations, but with the drawback of requiring 10 or 
more cameras. These algorithms are inspired by the geomet- 
ric observation that, from the optical center of each square - 
pixel camera, two lines can be identified in the projective 
reconstruction that must intersect the absolute conic. The 
absolute quadratic complex (AQC) encodes the set of lines 
intersecting this conic by means of a quadric in P^, which 
is the natural space containing Pliicker coordinates of lines. 
The AQC, being represented by a homogeneous symmetric 
6x6 matrix satisfying a linear constraint, depends linearly 
on 21 — 1 — 1 = 19 non-homogeneous parameters, which ex- 
plains the need for such a high number of 10 views. 

However, an informal parameter count reveals that far 
fewer cameras are theoretically sufficient. In fact, our main 
unknown is a space homography, which depends on 15 pa- 
rameters. Being our target reconstruction defined up to a 
similarity, which depends on 7 parameters, there are 15 — 
7 = 8 unknowns left to be found, which correspond to the 
degrees of freedom (dof) needed to determine the plane at 
infinity (3 dof) and the absolute conic within it (5 dof). Know- 
ing camera skew and aspect ratio amounts to two equations 
per camera and thus at least 4 cameras should be given in 
order to solve the problem. Given the non-linear nature of 
these equations, multiple solutions can be expected and so 5 
cameras should be the minimum required to obtain, generi- 
cally, a unique solution (see ['6, Table 19.3]). 

The main goal of this paper is to obtain a Euclidean re- 
construction from the minimum number of cameras using 



exclusively the pixel shape restriction. The geometric ob- 
ject that will be employed for this purpose is the variety 
of conies intersecting six given spatial lines simultaneously, 
which will be termed the six-lines conic variety (SLCV). 
The SLCV as a geometrical entity has been studied in ITSl . 
although our treatment is independent and self-contained. In 
this paper we are interested in the SLCV given by the abso- 
lute conic at infinity. 

The SLCV for six lines in generic position can be iden- 
tified with a surface of P^* (i.e., the projective space given 
by the planes of P^) of degree 8. We prove that this degree 
reduces to 5 in the case of the three pairs of isotropic fines 
of three finite square-pixel cameras. We show that the fifth- 
degree SLCV has three singularities of multiplicity three, 
given by the three principal planes of the cameras. This re- 
sult is used in Algorithm 1 (see Table |2|i to generate a bidi- 
mensional parameterization of the candidate planes at in- 
finity compatible with three square-pixel cameras. This pa- 
rameterization, together with the additional data given by 
another two or more square-pixel cameras permits to iden- 
tify the true plane at infinity through a two-dimensional op- 
timization process, leading to Algorithm 2 (see Table |3]l. 
However, the technique could as well use other additional 
data such as some scene constraints. 

Experiments with real images for the autocalibration of 
scenes with 5 or more cameras with square pixels and oth- 
erwise varying parameters are provided, showing the good 
performance of the proposed technique compared to other 
autocalibration methods. In the absence of knowledge about 
the principal point of the cameras, the SLCV algorithm turns 
out to be the only feasible approach to solve the autocalibra- 
tion problem in the minimal case of 5 cameras up to the 
case of 9 cameras. For 10 or more cameras, the results are 
similar to those of the AQC algorithm. Table [1] summarizes 
the similarities and differences between the autocalibration 
approaches based on the AQC and the SLCV 

The paper is organized as follows: The basic background 
for the problem is briefly recalled in section |2] Section |3] 
presents the SLCV along with the basic algebraic geometry 
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tools required for its definition and analysis as well as our 
main theoretical results. The algorithms motivated by these 
results are presented in section |4] and the corresponding ex- 
periments are shown in section |5] Conclusions of the paper 
are found in section |6] An advance of some of the results of 
this paper appeared in the conference paper [3J. 



and then 



X^ + Y^ +Z^ = {XJ,Z){X,Y,Z)'^ 



±i 



2 Camera model, square pixel cameras and projective 
calibration 

We suppose that the cameras can be described using the pin- 
hole camera model 161, which is defined by the optical center 
C and the projection plane endowed with an affine coordi- 
nate system given by the pixel structure of the camera. 

The equations of the projection are of the form x ^ PX 
where P = K(R| -RC), X == {X,Y,Z,Ty are the homoge- 
neous Euclidean coordinates of a 3D point, x = {x,y,z)'^ are 
the homogeneous coordinates of its projection, R is a rota- 
tion matrix, C are the inhomogeneous Euclidean coordinates 
of the optical center, and K is the intrinsic parameter matrix, 
given by 

^fnix —f nix cot 9 "0^ 
K= I /mv/sin0 vq 
1 

where / is the focal length, nix and m, are the number of 
pixels per distance unit in image coordinates in the x and y 
directions, 9 is the skew angle and {uo,vo) is the principal 
point. 

If the camera aspect ratio T = niy/ntx and the skew angle 
9 are known, the affine coordinate transformation 

^TCOS0 0^ 

sin 
,0 1, 

of the image plane permits to assume that the intrinsic pa- 
rameter matrix has the form 



This equation can be expressed as 



1^0)1 = 0, r«i = o, 



'a Uq 

K = I a v() 

,0 1 



(1) 



which is the intrinsic parameter matrix of a square pixel 
camera, i.e., one for which T —\ and 9 — n/2. 

The back-projected lines of cyclic points at infinity I ^ 
(l,j,0)^ and i = (1,— /,0)^ are the isotropic lines of the 
camera. These lines intersect the absolute conic, for if X = 
(X , F, Z, 0) ^ is the intersection of one of these two lines with 
the plane at infinity, we have 

(l,±/,0)^ - PX = KR(X,F,Z)^, 



(2) 



so that 



{x,Y,z)' -r'k-'(i,±/,o)', 



which is the square pixel condition in terms of the image of 
the absolute conic (lAC) O) == (KK^)"'. 

We recall here |l6l that it is possible to obtain a pro- 
jective calibration from image point correspondences only. 
This means that, given a sufficient number of projected points 
x,y obtained with Nc > 2 cameras, we can obtain a set of 
matrices P,- and a set of point coordinates X^ such that x,y ~ 
P,Xj, where P, = P,H^' and Xj = HX,- for some non-singular 
4x4 mati'ix H. 

Euclidean calibration can be defined as the obtainment 
of a matrix H changing the projective coordinates of a given 
projective calibration to some Euclidean coordinate system. 
It is well-known that Euclidean caUbration up to a scale fac- 
tor is equivalent to the recovery of the absolute conic at in- 
finity f2oo. 



3 The six line conic variety 

Let us suppose we have a projective calibration of three 
square pixel cameras (or, as explained before, cameras with 
known pixel shape) and let P, be their camera matrices and 
C, their centers, i~ l,2,3.LetusdenotebyZ,,/; the isotropic 
lines of camera /. The plane at infinity tToo will intersect these 
lines in points of the absolute conic. Therefore the planes 7t 
candidates to be the plane at infinity are those intersecting 
the isotropic lines in points of a conic (see Figure[T]l. We are 
going to see that these planes are given by a 5th-degree alge- 
braic equation in the coordinates 71^ {ui,U2,u^,U4)~^ of the 
plane. In order to obtain this equation, some mathematical 
preliminaries are needed. 



3.1 The equation of the six-line conic variety 

We recall that lines of P^ are in one-to-one correspondence 
with non-null singular antisymmetric 4x4 matrices defined 
up to a non-zero scalar factor The correspondence is given 
by the mapping that assigns to the line / passing through 
points p,q the Pliicker matrix L — M(p,q) = pq^ — qp^ ||6l. 
There is an equivalent mapping attaching to the line deter- 
mined by planes a,j3 the matrix L* = M(a,j3). These two 
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Result 1 Given six lines with Pliicker matrices L,- and vec- 
tors ;r,ai, . . . ,34 G C , let us consider the polynomial 



Fig. 1 Illustration of the incidence relations between the isotropic lines 
of three cameras, the plane at infinity and the absolute conic. 



matrices are related by the transformation * : L = ('«;;) H> L* 
where 

/ m34 m42 m23\ 

— m34 mi4 OT31 

— OT42 —mi4 mi2 
\-OT23 -'«31 -mi2 / 

We also recall that the intersection of the line of Pliicker 
matrix L with the plane n = {ui, M2, M3, M4)^ is the point p = 
Ln, which is zero if and only if the line L is contained in K. 

If H is a coordinate change X' = HX, the line L is written 
in the new coordinate system as 



(3) 



l' = hlh^. 

The degree-two Veronese mapping V2 : P" 

n(n + 3)/2, is defined by 

V2{xi,...,X„+l)'^ = ({xiXj)i<jY . 

In particular, for n = 2, 3 we define: 

V2(xi,X2,X3)^ = (xi,XiX2,XiX3,X2,X2X3,X3)''", 
V2(xi,X2,X3,X4)^ = {x\,X\X2,XiXT,,XiXA,x\, 
X2X3 1 ■^2-^4 7 -^3 I X^X4 , X4 ) 

Observe that a point x G P^ belongs to the conic of ma- 
trix C — (cij) if and only if 

/ \Tr. n u r, f ^12 Ci3 C23 \T 

V2(x) C = 0, where C= (^cii,—,—,C22,—,C33J . 

Therefore the points qi, . . . ,q6 G P^ He on a conic if and 
only if 

det(v2(qi)---V2(q6))=0. 

Analogously, 10 points q, lie on a quadric of P-' if and only 
if 

det(v2(qi),...,V2(qio))=0. 

Next we use the Veronese mapping to characterize whether 
the intersection of lines with a plane is contained in points 
of a conic. 



D(;r,ai,...,a4) ^ 

det (v2(Li;r), . . . , V2(L6;r), V2(ai), . . . , V2(a4)) . 



(4) 



The set of planes K intersecting the six lines in points of a 
conic is defined by 

D(;r,ai,...,a4) =0/orfl//ai,...,a4 G C*. 

Proof First note that if D(;r, ai , . . . , a4) == 0, the points L,-;r, ay 
lie on a quadric Q. By construction, the points L,-;r are also 
on the plane n, and therefore they lie on the conic Kr\Q. 

Conversely, let us suppose that the points L,7i; lie on a 
conic C contained in n. Let us choose coordinates [x^y^z^t) 
such that the plane n is given by the equation t — Q. Let 
C(x,y,z) = f = be the equations of C. Any quadric Q con- 
taining C has an equation of the form 

aC{x,y,z)+t{ax + by + cz + dt) = 0. 

Given four points Aj = {xj,yj,Zj,tj)'^ , it is always possible 
to find some quadric Q through them because the points lead 
to a hnear homogeneous system of four equations in the five 
unknowns a,b,c,d, a which always admits at least one non- 
trivial solution. Since the 10 points L,7r,ay lie on Q, we con- 
clude that D(;r,ai,...,a4) = 0. D 

Being each column V2(L,7r) of degree two in n, the equa- 
tion D(;r,ai , . . . ,a4) = is of degree 12 in the coordinates 
of the plane. Next result shows that we can factor out four 
trivial linear factors and obtain an 8th-degree polynomial in 
7t which is also shown to be independent of the Sj. 

Result 2 The set of planes % intersecting the six lines with 
Pliicker matrices L,- in points of a conic is given by the 8th- 
degree polynomial equation 



F{n) = Q 

defined by the relationship 

D(;r,ai,...,a4) = 

det(ai, . . . ,a4) (;r"^ai) ■ • • {n^aA)F{n), 



(5) 



(6) 



Furthermore, the polynomial F does not depend on the vari- 
ables aj. 

Proof Any plane 7t through the point ai is a trivial solution 
of D = 0, since the seven points LiTT,. . . ,Lg7r,ai would lie 
on the plane n and therefore the ten points lie on the quadric 
given by n and the plane formed by points a2 , a3 , a4. Conse- 
quently, the planes through ai produce a linear factor a:^ai 
of D and so do the planes through each of the three remain- 
ing points aj. 

Let us denote A = (ai, . . . ,a4). Next we show that det A 
also divides D. Since the determinant det A is irreducible |[T] 
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p. 176] when regarded as a polynomial in the coordinates of 
the points ay, it is enough to show that D vanishes when- 
ever det A = 0, i.e., when the points a^ are contained in some 
plane n' . In such case, the ten points L,7r, i ~ 1, ... ,6 and 
aj, y = 1 , . . . ,4, lie in the degenerate quadric ;r ■ ;r' so that D 
cancels. We have therefore the required factorization (0. 

Let us now show that F does not depend on the vari- 
ables Sj. The degree of the variables aj on the left hand side 
of ^ is two because the vectors V2(aj) are homogeneous 
quadratic in the entries of ay and each of the V2(ay) appears 
only once as a column of the determinant. On the right hand 
side of (|6|l we have the factor det A, which is homogeneous of 
degree one in each Sj, and the factors ?r^ay, which are also 
homogeneous of degree one. Therefore the remaining factor 
F does not depend on the a/ since otherwise there would be 
a mismatch between the degrees of both sides of (|6]l. 

Finally, let us check that the equation F{k) =0 char- 
acterizes the planes intersecting the six lines in points of a 
conic. Now that we have proven factorization (|6), it is trivial 
thatifF(;r) =0thenD(a:,ai,...,a4) = and, by Result[T] ;r 
is a plane intersecting the six lines in points of a conic. Con- 
versely, if n intersects the six lines in points of a conic, the 
determinant D vanishes for any ay. In particular, choosing 
any non-coplanar points ay ^ 7i, so that detA 7^ ^ Tf^a/, 
we conclude that F{n) =0. D 

The surface of P^ given by the planes intersecting the 
six lines in points of a conic will be called the six-line conic 
variety (SLCV). 



3.2 The SLCV given by the isotropic lines of three finite 
square-pixel cameras 

In the particular case of three finite square-pixel cameras, 
the configuration of the isotropic lines, which intersect pair- 
wise in the optical centers C,-, allows to further simplify the 
factorization from Result |2] 

Result 3 Given a projective reconstruction for three finite 
cameras, the planes at infinity compatible with the square- 
pixel property of the cameras are given by the real zeros of 
a 5th-degree polynomial G given by the following equation: 



F{k) = {n^Cx){n^ C2){n'^Ci) G{n). 



(7) 



Proof We have proved in Result |2] that the set of planes in- 
tersecting the lines L, in points of a conic is given by the 
zeroes of an 8th-degree polynomial, F{n). In the case of six 
isotropic lines of three square-pixel cameras, there are triv- 
ial solutions of the equation F{n) = which are of no inter- 
est, namely those corresponding to planes passing through 
any of the optical centers C,, i — 1,2,3. In fact, any plane 
through any of the C, intersects the six lines in at most 
five different points. Since five points always lie on a conic. 



all the planes of the three stars of planes through the opti- 
cal centers are zeros of F. Moreover, these are not candi- 
date planes at infinity since we are assuming finite cameras. 
Therefore, we can further factorize the polynomial F as in 
equation (|7]i. D 

In the context of camera autocalibration, we will also 
use the term SLCV for the 5th-degree variety, from which 
the three trivial linear factors have been removed. 

A straightforward application of this result is the direct 
obtainment of the candidate planes at infinity, when two 
vanishing points are known, as the real solutions of a fifth- 
degree polynomial equation in one variable (see [3]). 

Corollary 1 If two points of the plane at infinity are known, 
there are at most five candidate planes at infinity which can 
be found solving the 5th-degree equation in the homoge- 
neous coordinates A ; jX 

where a and j3 are any two planes intersecting in the line 
determined by the two points at infinity. 



3.3 Singularity of the principal planes 

From now on, we suppose as above that the six lines are the 
isotropic lines of three square-pixel cameras. The following 
result is the key to find a computationally efficient way to 
parameterize the set of associated candidate planes. For our 
configuration of three cameras, we will say that a principal 
plane is generic if it does not contain the projection centers 
of any of the other two cameras. 

Result 4 Any generic principal plane is a singularity of mul- 
tiplicity three of the variety of candidate planes G{k) — 0. 

Result m suggests to parameterize the set of candidate 
planes considering the pencils of planes with base contained 
in one of the principal planes. We will denote by Ujt the 
pair of lines through the optical center C,- and L,, L, their 
Pliicker matrices. The corresponding principal planes will 
be denoted as Kt, i = 1,2,3. Let us assume that n\ is a 
generic principal plane. To each pencil of planes that in- 
cludes n\ we will associate two candidate planes. To param- 
eterize the pencil, we can consider its element ^ through C2, 
so the pencil is given by the planes of the form Xk\ + Ii£, 
(see Figure|2|l. 

Result|4]guarantees that the 5th-degree polynomial G{Xni 
jiE,) factorizes as 

G(A;ri+jti<^)-M^//(A,M) (8) 

where // is a homogeneous polynomial of degree 2 whose 
zeroes provide the two candidate planes associated to £, . The 
approach of using a triple point to parameterize a 5-th degree 
variety is illustrated in Figure |3]for the case of a 2D curve. 
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Fig. 2 Elements in the parameterization of the SLCV given by Re- 
sultE] 



Proof [Proof of Result m Using Result|2]we have that 

D = det (v2(Li7i;), . . . , V2(L37r), V2(ai), . . . , V2(a4)) 
= det(ai, . . . ,a4) {%^ai) ■ ■ ■ (k'^m) (;r"^Ci)(;r"^C2) 
(7i;^C3)G(7r). 

With the notation introduced above, letn — Xni+iJ.£,. Since 
hji are contained in tti we have LiTTi = LiTTi =0 and 
therefore LiTT = /xLi^ and Lin = jiLiS,. Consequently 

D(A,Ai)=M^det(v2(Li^),V2(Li^),V2(L27r),...,V2(a4)). 

On the other hand, since Ci G Ki, we have Cjn = jxCj ^, 
so 

;u'*det(v2(Li^),V2(Li^),V2(L2;r),...,V2(a4)) 
- Aidet(ai , . . . ,a4) (71^31) • • • (n'^a^) (^^Ci)(7rTC2) 

{7t'^C3)G{7t). 

Using the genericity of tti and choosing conveniently the 
points ay we have that C2,Cj,,aj ^ Ki for j — 1, ... ,4, i.e., 
jX does not divide any linear factor besides Tt^C\. Therefore 
IV" divides G{X%\ + ^it,) and so Tti is a singular point of G 
of multiplicity three. 

D 



4 Algorithms 

We have seen that the set of planes compatible with three 
cameras with square pixels is a surface given by a fifth- 
degree equation 

This will allow to compute the plane at infinity by means 
of a two-dimensional search if enough additional data about 




Fig. 3 Illustration of the approach to obtain the points of the SLCV. 
The projective plane curve -F(p) =0 is of degree five and has po as 
a triple point. Therefore, substituting in this equation the parametric 
equation of a line through this point, results in an equation of the form 
F(Jlp() + /iqo) = fl^G{X,fl), so that the other two points of intersection 
of the line and the curve are obtained by solving the quadratic equation 
G{A,;U) = 0. 

The SLCV is a fifth-degree surface in dual space P^* and the proposed 
technique to parameterize is based in computing its intersections with 
lines through one of its triple points. These lines in P^^* correspond to 
pencils of planes in P^^ . 



the cameras or the scene are available. In the first part of this 
section we will address a convenient way to parameterize 
the surface of candidate planes, while in the second part we 
will focus on the particular case in which the additional in- 
formation stems from the presence of two or more auxiliary 
square-pixel cameras. 



4. 1 Parameterization of the candidate planes 

Exploiting Result|4]and assuming that the first principal plane 
Tti is generic, our algorithm will sweep the set of candidate 
planes at infinity corresponding to three square-pixel cam- 
eras. We will define a one-to-two mapping attaching to each 
real line / of ttj the two candidate planes Xi, X2 containing 
/, i.e., the intersection of the pencil of planes through I with 
the variety of candidate planes. Hence, if we parameterize 
the set of real lines / of ;ri that do not contain the optical 
center Ci, we will obtain accordingly a two-fold parame- 
terization of the set of candidate planes. To this purpose we 
fix a point r in /i which, together with the optical center Ci 
will parameterize the points q G /i as q = r + zCi, being z 
a complex numbeo We now define / as the line through q 
and q (see Figure|4]l. 

We will consider as generators of the pencil of planes 
through / the principal plane K 1 and a plane ^ passing through 

' The fixed point r may be chosen as the intersection of /] with the 
plane whose coordinates coincide with those of the point C 1 . 
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I and C2, so the planes of the pencil are of the form 

As explained in the previous section, the solutions con- 
tained in the pencil through / will be the zeroes of the poly- 
nomial H{X,li) defined in ([8]). Since it is a homogeneous 
degree two polynomial, it has an expression of the form 



H{X,l^)=AX^+BXfi 



Cpl\ 



(9) 



where A^B, C are polynomials in the coefficients of the Pliicker 
matrices L,, L/. Their size, of the order of thousands of terms, 
is not suitable for algorithmic use. However, a convenient 
coordinate change shortens them so that they add up to about 
one hundred terms. 

Specifically, denoting by q and q the intersection points 
of the isotropic lines l\ and li with the plane ^ and assum- 
ing that the points C,, q and q are in general position, we 
consider the coordinate change H performing the mapping 

Ci^(0,0,0,l)T 
C2^(0,0,1,1)T 

C3^(0,1,-1,1)T (10) 

q>^(l,/,0,0)T 
q^ (1,-^,0,0)^. 
In these new coordinates we obtain a polynomial 



//o(A,M) =AoA2 + BoAm + CoM^ 



(11) 



where the coefficients Aq , Bo i Q depend only on L2 , L3 , since 
Li is constant in the new coordinate system. 

In Algorithm 1, we assume that the expressions of Aq, Bq 
and Co have been precomputed (see the appendix for further 
details). 



4.2 Euclidean calibration with five or more cameras 

The parameterization provided by the previous Algorithm 1 
can be employed, in particular, to perform a two-dimensional 
search for the plane at infinity using the knowledge provided 
by two or more additional square-pixel cameras. The next al- 
gorithm explores the set of solutions associated to the first 
three cameras, aiming at the minimization of a cost function 



C(z)-min{CoUi(z)),Co(Z2(z))}, 



(12) 



where C(i{x) measures the compatibility of candidate plane 
X with the square pixel condition of the additional cameras. 
For a given x^ the lAC of the cameras, C0i{x), is calcu- 
lated as the conic through the projections of the points of in- 
tersection of Z, andZ^, i~ 1,2,3, with X- Then the cost Co (x) 

^ See appendix |B] for explicit formulas of Pliicker matrices of the 
isotropic lines. 



Algorithm 1 

Objective 
Given: 

1. a projective calibration of three cameras with optical centers C; 
and isotropic lines /;,/) with Pliicker matrices L;,L;, with generic 
principal plane Tti , 

2. a fixed point r 7^ Ci on the isotropic line l\ and 

3. a complex number z, 

compute the two planes through the points q = r + zCi 6/1 and 
q = f + zCi 6 f 1 that intersect isotropic lines hjijhji in points 
that, together with q and q, lie on a conic. 
Algorithm 

1. With q = r + zCi compute the matrix of the coordinate 
change jlOt and apply ^ to transform L2 and L3 . 

2. Using precomputed expressions, calculate the coefficients of 
the second-degree homogeneous polynomial l IlU and its roots 
(A,-,AJ,),i=l,2. 

3. Obtain the solution planes in the transformed coordinates as 
X: = A,?ri +Ai;^, with Xi = (0,0, 1,0)^, ^ = (0,0,1,-1)^, 
( = 1,2 and revert coordinate change for the planes Xi' thus 
obtaining solution planes Xiiz), Xii^)- 



Table 2 Algorithm to compute two planes ;i;i(z), Xii^) intersecting 
the isotropic lines of three square-pixel cameras in points of a conic, 
where z parameterizes the set of lines in Tti. 



is computed from these lACs. In the course of the algorithm, 
complex solutions may arise, albeit the actual lAC must be 
real. Therefore, some additional constraints are taken into 
account in the design of Co{x)- In particular, the cost Co(;i;) 
is the maximum of the weighted sum of four non-negative 
terms: 



Coix) = max ^ 7*Q.(«,(z)), 



(13) 



where the weights Yk > 0, Ci (©) penalizes complex solu- 
tions, C2{co) discourages non positive-definite lACs, Cj,{co) 
measures the square pixel condition (l2) and C4{co) penalizes 
principal points outside the image domain. 

Before computing the individual costs, the lACs co un- 
dergo two normalization steps. First, the homogeneous ma- 
trix CO is scaled by the unit complex number s that maxi- 
mizes the Frobenius norm of 9{{ia)}. This is a constrained 
optimization problem whose solution is given by a biquadratic 
equation in 3is. Then, co is scaled to unit Frobenius norm, 

II-IIf. 

Next, let us describe each term in (fT3T l. If u = vec{3ico) 
and V — vec(3 co) are the vectorizations of the upper- triangular 
part of the real and imaginary parts of CO, 

C,(a))^||uvT-vuT||f/(||u||2-K||v||2). 

The term C2 is motivated by Sylvester's criterion, which 
states that a hermitian matrix is positive-definite if and only 
if all of its leading principal minors are positive. Let D, (A) , / = 
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Fig. 4 Geometric elements used in Algorithm 1. The three square-pixel cameras are represented by their optical centers Cy and their isotropic 
lines Ij, tj contained in the principal planes Ttj, j = 1,2,3. Planes Tti and ^ are generators of the pencil of planes through the line joining q(z) and 
q(z)- Within the pencil, solution planes Xiiz), i= 1,2 intersect the six isotropic lines in points of a conic. 



1 , 2, 3 be the three leading principal minors of matrix A and 
^(A) = -lf=imin{0,A(A)},then 

C2(«)=minU(9^{a)}),g(9^{-«})}, 

penalizes the non positive-definiteness of the homogeneous 
matrix o. 

The term C3 measures the deviation from the square pixel 
condition. We choose: 



C,{(0\ 



^JT^-ll+COS^e, 



where T — niyjin^ = (Oi 1/(022 and cos^ 9 — 0)^2/(^11^2)- 
Observe that Ci{cOi) = for / = 1,2,3, by construction of 
the SLCV, but not necessarily for the additional cameras. 
The two terms in C3 may be weighted differently, too. 

The term C4{co) is the taxicab distance (|| • ||i) from the 
principal point {uo,vo) to the boundary of the image domain 
if the principal point lies outside the image domain and zero 
otherwise. Formulas for uo,vo in terms of (0 are well known 

70)33, where 



(0 



'23/ "'SS 



in the literature: mq = 0)^*3/0)33 and vq 

CO* ^ 0)^'. This term accounts for spurious solutions with 

unrealistic location of the principal points of the cameras. 

Algorithm 2 performs a nonlinear optimization of cost 
function (fT2l i in two steps. In the first one, the complex plane 
is sampled and the sample of minimum cost is selected. In 
the second one, a search for a local minimum is performed 
using as starting point the complex value provided by the 
first step. 

A sampling of the complex plane that has been empiri- 
cally found to be useful consists in splitting the plane into 
the unit disk and its complement, and employing a uniform 
sampling in modulus and phase for the first, and the inverses 
of these values for the second. The Nelder-Mead method 
(downhill simplex method) outperforms other optimization 



algorithms with numerical first derivatives such as steepest 
descent or conjugate gradient. 

The Euclidean upgrading may be computed from the 
stratified approach in Algorithm 10.1 of [6|, with the plane 
at infinity and the lAC of one of the cameras provided by 
Algorithm 2. 



5 Experiments 

5.1 LED bar dataset 

The proposed technique has been tested in the calibration of 
a set of five synchronized square-pixel video cameras with a 
resolution of 1280 x 960 pixels using a bright point device. 
Instead of a single-point device, a rigid bar with three light- 
emitting diodes (LEDs) is employed in the tests. In addition 
to providing ground truth (the bar length constancy) to test 
the quality of the results, this allows to use for comparison a 
calibration algorithm based on the geometry of the captured 
set of 3D points. 

In the experiments, a projective reconstruction of the 
scene is first obtained for the five cameras. This is accom- 
plished using Algorithm 10.1 of [6] (eight-point algorithm) 
to compute the fundamental matrix of a pair of cameras 
and then alternating 12.2 (triangulation) and 7.1 (resection) 
of 16J. The result is then optimized by projective bundle ad- 
justment ||6l p. 434]. 

Euclidean upgrading is performed using three different 
techniques: Algorithm 2, an algorithm based on the dual ab- 
solute quadric (DAQ), and an algorithm based on the LED 
bar geometry. 

The Euclidean upgrading algorithm based on the DAQ 
assumes square pixels and principal point at the center of 
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Algorithm 2 

Objective 

Given a projective calibration of W^. > 5 cameras with square pixels, 

projection matrices P, and isotropic lines /,,/;, obtain a Euclidean 

upgrading. 

Algorithm 



in 

.,N- 



the 
= 1,...,MJ 

l,k=\,.. 



,M>, compute cost 



1. For each z in the set Z = {0} 

C{z) as follows. 

(a) Compute planes X\{z) and Xii^) using Algorithm 1. 

(b) For each plane ;tG{Zi (2)^2(2)}. 
i. Compute the points of intersection of isotropic lines 

/;, li,i= 1,2,3 with the plane, project them with pro- 
jection matrix Pi and obtain the matrix ft)] of the 
conic that they define. This is the lAC of the first 
camera. 

ii. Obtain the lACs (B, for cameras i = 2,...,N^. by 
transferring ff)i onto them using the plane X'^i ^nd 
P; (Note that the same result would be obtained if 
the second or third camera was employed). Normal- 
ize them to (unit Frobenius norm and maximum real 
part). 

iii. Compute Co{x), according to equation J13t . 

(c) Compute C(z) =min{Co(;i;,(z)), Co(;i;2W)}- 
Take zo = arg min^gz C(z) . 

Perform non-linear optimization to obtain a local minimum z\ 
near the value zo, repeating steps (a), (b) and (c) above for each 
evaluation of the function. Choose, of the two planes attached 
to zi, the one with minimum cost Cq. The Euclidean upgrading 
is determined by taking this plane as the plane at infinity and 
the associated conic as the absolute conic. 



Table 3 Algorithm to compute a Euclidean upgrading from a projec- 
tive calibration. 



After the Euclidean calibration, a Euclidean bundle ad- 
justment is performed, including the enforcement of the square 
pixel shape and taking into account the lens distortion co- 
efficients (according to the four-parameter OpenCV model 
0). 

The sampling parameters in Algorithm 2 have been M = 
A^ = 50, so that the cost function has been evaluated on 2 x 
50 X 50 = 5000 points. Figure|5]shows a sample of the input 
data for the calibration process. Figure |6]shows two views of 
the sampled points of the SLCV computed in Algorithm 2. 
Figure|2]displays the cost function (fT2l l for this experiment, 
normalized and pseudo-colored from blue (0) to red (1). 
A logarithmic transformation has been applied to enhance 
the visualization of small costs. The reconstructed scene ob- 
tained with Algorithm 2 is shown in Figure |8] Table|4]shows 
the reprojection error and the quotient between the standard 
deviation and the average of the segment lengths for each of 
the Euclidean upgrading techniques, before and after bundle 
adjustment. 

Input data (triplets) - camera 




200 400 600 800 1000 1200 



Input data (triplets) - camera 2 



the image. This provides four linear constraints on the DAQ 
for each camera f2T\. 

The algorithm using the LED bars is based on the fact 
that three aligned equidistant points determine the point at 
infinity of the common line ||6| p. 50]. Thus each captured 
position of the LED bar provides a point at infinity in the 
projective reconstruction, and, therefore, three positions of 
the bar in the projective reconstruction determine the plane 
at infinity and, consequently, an affine calibration. Euclidean 
calibration is then obtained by computing the affinity that 
makes the lengths of the segments coiTesponding to cap- 
tured positions of the ligid bar as close as possible to the 
true bar length. This requires at least six captures of the bar 
and the solution of a least-squares problem followed by a 
Cholesky factorization f20l to determine the affinity up to a 
Euclidean motion. Affine and Euclidean calibration are re- 
peated 500 times using different random sets of three bars 
for the affine calibration and all the bars for the Euclidean 
calibration, and the reconstruction with smallest bar length 
variance is selected. 




200 400 600 800 1000 1200 



Fig. 5 Example of input data for the tests (triplets of aligned LEDs in a 
rigid bar). Only the points are used in the calibration with the proposed 
algorithms. The triplet structure is used for evaluation puiposes. 



5.2 Checkerboard dataset 

The developed method has also been tested on a different set 
of 5 to 10 images that contain checkerboard patterns [114l . 
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DAQ 


SLCV 


Triplets | 


Euc. calib. 


Euc. BA 


Euc. calib. 


Euc. BA 


Euc. calib. 


Euc. BA 


Rep. error 


15.04 


1.74 


1.71 


0.64 


13.56 


0.63 


ff/M 


0.15 


0.029 


0.049 


0.0098 


0.026 


0.0048 



Table 4 Results of experiment comparing calibration techniques. The common previous projective calibration shows a reprojection error of 1 .23 
pixels. <y and fl are, respectively, the typical deviation and the average of LED bars lengths. 




0.3 
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0.1 
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-0.2 
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Fig. 6 Two views of the sampled points of the SLCV. 



Fig. 8 Two views of the 3D reconstmction in one of the experiments. 




-1 -0.5 




I 






Fig. 7 Experiment with the LED bar dataset: plots of the sampled cost 
function. Left: values at the unit disc of the complex plane, |z| < 1. 
Right: values at the complement of the unit disc, at positions 1/z. The 
white cross in the left plot marks the location of the minimum. 



This calibration rig provides an additional validation of the 
results. The images, of size 1280 x 960 pixels, were acquired 
with a Sony DSC-F828 digital camera. To test varying pa- 
rameters, the equivalent focal length (in a 35 mm film) of 
the camera was set to 50 mm in eight of the images and to 
100 mm in the remaining two. Variations due to auto-focus 
were not controlled. For this range of focal lengths, the lens 



distortion (radial and tangential) can be neglected, so it is 
assumed to be zero. 

For each of the following experiments, scale-invariant 
key points (SIFT |[8|) are detected and matched across im- 
ages, obtained using ifTTl . Then, a projective reconstruction 
of the scene is obtained, as already explained. Normalization 
of coordinates ("preconditioning") is applied, as it is essen- 
tial to improve the numerical conditioning of the equations 
in the different estimation problems involved (fundamental 
matrix, triangulation, resection and bundle adjustment). Ta- 
ble 15.21 summarizes the parameters of the projective recon- 
structions. The resulting set of projection matrices are the 
input to the autocalibration algorithms. 

Table|6]shows the results for the experiment with 10 im- 
ages. In this case, it is possible to compare the SLCV and 
DAQ algorithms to a technique based on the AQC ll22l . A 
Euclidean bundle adjustment with enforcement of the square- 
pixel condition is performed after the Euclidean calibration 
with each of the four compared techniques. The solution ob- 
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Method 


AQC 


DAQ 


SLCV 


Euc. BA 


Mean focal length a 


1822.9 
3433.2 


1851.4 
3579.3 


1806.2 
3399.4 


1847.8 
3555.7 


a standard deviation 


10.95 
36.15 


4.86 
13.33 


16.1 

29.21 


7.77 
4.53 


Meanpp (mq.I'o) 


(623.1,463.9) 
(618.6, 446.8) 


(638.1,477.8) 
(654.5, 505.7) 


(607.4, 468.9) 
(605.5, 455.4) 


(620.4, 485.3) 
(588.0, 525.2) 


p.p. std. dev. 


(16.85,8.13) 
(34.44, 0.55) 


(1.61,4.96) 
(17.27, 10.91) 


(27.7, 11.15) 
(53.18,7.47) 


(6.72, 8.71) 
(24.79, 6.36) 


Rep. error 


0.65765 


0.27181 


0.27497 


0.22386 



Table 6 Reprojection eiTor and intrinsic parameter comparison for the experiment with 10 images of the Checkerboard scene. For each statistic, 
the top row corresponds to the value for cameras with / = 50 mm (equivalent in 35 mm film) and the bottom row corresponds to cameras with 
/ = 100 mm. Data are given in pixels. 



Scene 


Checkerboard 


Plaza de la Villa | 


# images 


10 


5 


16 


5 


# 3-D points 


1494 


1199 


8555 


1083 


# image points 


6660 


4472 


35181 


3483 


Rep. error (BA) 


0.22268 


0.17678 


0.16053 


0.17855 



Table 5 Projective reconstructions. 





Fig. 9 Experiment with 10 images of the Checkerboard dataset: plots 
of the sampled cost function. Left: values at the unit disc of the com- 
plex plane, \z\ < 1. Right: values at the complement of the unit disc, at 
positions 1 /z. The white cross in the left plot marks the location of the 
minimum. 



tained (regardless of the AQC, DAQ or SLCV initialization), 
is given in the last column of table |6] 

If the image resolution and the size of the CCD of the 
camera are known, it is possible to convert the focal length 
from mm to pixels. For this dataset, an equivalent focal length 
of / = 50 mm translates to a = 1850 pixels, which is very 
close to the values obtained by the different algorithms tested 
(first row of table |6]l. 

Because the DAQ algorithm yields good estimates of the 
intrinsic parameters (with small dispersion and reprojection 
error), we may conclude that the principal point (p.p.) of the 
camera is close to the center of the image. This observation 
is also supported by the estimates of the p.p. due to the other 
algorithms. All three autocalibration methods show a strong 
agreement with the three-homography calibration algorithm 
inL14JandL6,p.211]. 

Euclidean calibration by means of the developed tech- 
nique (SLCV) provides competitive results with respect to 
the other methods (DAQ or AQC), but with a slightly bigger 
dispersion around the mean values. In normalized coordi- 
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Fig. 10 Experiment with 10 images of the Checkerboard dataset: re- 
constructed 3D scene. 



nates, the weights used in ( flj] ) to measure the goodness of 
fit of the lACs in Algorithm 2 are 71 = 72 = 73 = 74 = 1 ■ The 
sampling parameters in Algorithm 2 are M = A^ = 100. Fig- 
ure 15.21 shows the sampled cost function. Figure [TO] shows 
two views of the densely reconstructed scene correspond- 
ing to the Euclidean calibration by means of the SLCV al- 
gorithm. The dense reconstruction was obtained by feeding 
the images and the Euclidean calibration of the cameras to 
the Patch-based Multi-view Stereo Software (PMVS) soft- 
ware fsl- 

Next, an experiment with a subset of 5 images is car- 
ried out. To account for varying parameters, three of the im- 
ages correspond to a focal length / = 50 mm and the re- 
maining two have / = 100 mm. Table 15.21 shows the pa- 
rameters of the projective reconstruction. Table ^compares 
the Euclidean upgrading given by the DAQ and the SLCV 
algorithms. Both initializations converge to the same solu- 
tion after Euclidean bundle adjustment (last column of Ta- 
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Method 


DAQ 


SLCV 


Euc. BA 


Mean focal length a 


1833.3 
3511.2 


1855.2 
3602.1 


1883.1 
3662.9 


a standard deviation 


0.46 
16.24 


3.84 
7.90 


4.78 
26.15 


Mean p.p. (mo,vo) 


(637.4, 478.9) 
(664.2, 493.3) 


(624.5, 477.2) 
(618.0, 492.4) 


(620.0, 500.2) 
(588.8, 569.5) 


p.p. std. dev. 


(0.72,2.12) 
(28.61, 15.91) 


(1.41, 1.48) 
(31.53,20.03) 


(8.72, 0.68) 
(12.3, 15.4) 


Rep. error 


0.26204 


0.20134 


0.17689 



Table 7 Reprojection error and intrinsic parameter comparison for experiment with 5 images of the Checkerboard dataset. Same notation as in 
Table|6] Both DAQ and SLCV methods yield good results. 




^ 
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Fig. 11 Experiment with 5 images of the Checkerboard dataset: reconstructed 3D scene. 



ble|7J. The SLCV method yields similar results to those of 
the DAQ method but requiring fewer equations per camera. 
The reconstructed scene corresponding to the Euclidean cal- 
ibration by means of the SLCV algorithm is shown in Fig- 
ure [TT] The SLCV method provides a sensible initiaUzation 
to Euclidean bundle adjustment. 



Next, we demonstrate the good performance of the de- 
veloped method in case of decentered principal point. To 
do so, the previous 5 images of the Checkerboard dataset 
are extended to 1600 x 1200 pixels from the upper-left cor- 
ner. The principal point is, as seen in Table |6] near the point 
with coordinates (640,480) pixels, which significantly dif- 
fers from the new image center at (800,600) pixels. A pro- 
jective reconstruction of the scene is obtained, with 1197 3- 
D points, 4463 image projections and a reprojection error of 
0.17699 pixels. Table |8] compares the Euclidean upgrading 
by the two autocalibration methods in Table Q Because the 
hypothesis of known principal point (e.g. at the image cen- 
ter) is not satisfied, the DAQ method performs poorly. The 
SLCV method, however, yields the similar good results as 
in Table |2] because it solely relies on the square-pixel con- 
straint. The last column of Table |8] shows the result after 
refining the SLCV Euclidean calibration by bundle adjust- 
ment. 



5.3 Outdoor dataset 

As an additional experiment, 16 images of the Plaza de la 
Villa in Madrid were acquired with an Olympus E-620 dig- 
ital camera at a resolution of 1280 x 960 pixels. The focal 
length was set to / = 50 mm ia = 1850 pixels) in half of 
the images and to / = 70 mm (a = 2590 pixels) in the other 
half. Reconstructions were carried out with all images and 
with a subset of 5 images (with different focal lengths). The 
parameters of the projective reconstructions are summarized 
in Table 15.21 Figures fT2] and [T3] show reconstructions of the 
scene with 16 and 5 images, respectively, obtained by Eu- 
clidean upgrading with the SLCV method. 



6 Conclusions 

In this paper we have proposed an algorithm to obtain a Eu- 
clidean reconstruction from the minimum possible number 
of cameras with known pixel shape but otherwise varying 
parameters, i.e., five cameras. To this purpose we have intro- 
duced, as our main tool, the geometric object given by the 
variety of conies intersecting six given spatial lines simul- 
taneously (the six-line conic variety, SLCV). We have pre- 
sented an independent and self-contained treatment includ- 
ing a procedure for the explicit computation of the equation 
of the SLCV defined by three cameras. 
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Method 


DAQ 


SLCV 


Euc. BA 


Mean focal length a 


3143.7 
25499.9 


1827.2 
3516.0 


1877.6 
3643.56 


a standard deviation 


29.28 
20080.2 


3.87 
5.30 


5.37 
22.69 


Mean p.p. {uqvo) 


(801.3,516.4) 
(12192.1,-3639.6) 


(622.6, 478.5) 
(626.0, 493.9) 


(617.5, 500.6) 
(584.0,571.8) 


p.p. std. dev. 


(261.1,44.54) 
(23167.2, 9700.0) 


(6.49, 0.92) 
(43.46, 7.4) 


(7.54, 0.54) 
(15.97, 9.37) 


Rep. eiTor 


3324.95 


0.20725 


0.17707 



Table 8 Reprojection error and intrinsic paraineter comparison for experiment with 5 expanded images of the Checkerboard dataset. Same notation 
as in Table [6] The SLCV method clearly outperforms the DAQ method if the principal point is not near the image center 
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Fig. 13 Experiment with 5 images of the outdoor dataset: reconstructed 3D scene. 



While the SLCV for six lines in generic position is, in 
general, a surface of degree 8, we have shown that this de- 
gree can be reduced to 5 in the case of the three pairs of 
isotropic lines of three finite square-pixel cameras. We have 
seen that a direct application of this result is the obtainment 
of the candidate planes at infinity, when two points at infin- 
ity are known. 



We have seen that the fifth-degree SLCV has three sin- 
gularities of multiplicity three. We have used this fact to 
obtain a computationally efficient parameterization of the 
SLCV that, if we have some additional data, permits the ob- 
tainment of the plane at infinity by means of a two-dimensional 
search. An algorithm has been proposed for the case in which 
two or more additional cameras are available. 
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Fig. 12 Experiment with 16 images of the outdoor dataset: recon- 
stiTtcted 3D scene. 



Experiments with real images for the autocalibration of 
scenes with 5, 10 and 16 cameras have been given, show- 
ing the good performance of the SLCV technique. Particu- 
larly, in the case of fewer than 10 cameras, so that the AQC 
cannot be used, the SLCV method overcomes the limita- 
tion of the DAQ based technique for camera autocalibration 
with varying parameters, i.e., the need for known principal 
point. Thus, the developed method seems to be the only one 
that solves the autocalibration problem in the above situa- 
tion without requiring previous initialization. Furthermore, 
we have shown that the SLCV method is also a sound alter- 
native to other approaches that require 10 or more cameras 
or some knowledge of the principal point. 



Pliicker matrices of the lines /, are 
/OOO -1\ 
000 - 



Li=Ciq^-qC/ 



L2 = C2P2 — P2C2 -- 



L3=C3P3"-P3CJ = 



000 

V 1 / / 

/ -X2 -X2 \ 

->-2 -yi 

X2 yi -Z2 

\x2 yi Z2 J 

I -X3 X3 -x-i\ 

Xi Zi+yi ->'3 

-X3 ->'3 -Zi -Zi 
\ Xi J3 Z3 / 



where the points p, = {xi,yi,Zi,0)^ are the intersection of lines /,- with 
the plane K^ = (0,0,0, 1)^. The Pliicker matrices L, are just the com- 
plex conjugate of the matrices L,. Substituting in equation (|4} and fac- 
toring out the trivial linear factors we obtain the polynomial Ho{a,P). 
Explicit formulas can be found in 
http://www.gti.ssr.upin.es/~jir/SLCV.html. 



B Calculation of the isotropic lines 

The isotropic lines are the back-projection of the cyclic points at infin- 
ity (1,±!,0)^. Let us see how we can compute their Pliicker matrices 
from the rows of the con'esponding projection matrix P. In a projective 
reference in which P = (pi,p2,P3)^, the back-projection of the cyclic 
points is given by those X 6 P-' such that 

px~(i ±;o)^. 

Therefore the cross-product 

pXx (1 ±/0)'^=0, 

or, equivalently, X satisfies the equations 

P3TX = 0=(pJ±/p7)X. 

Hence, the isotropic lines are defined by the intersection of the planes 
P3 and p2 ± ;pi ■ Finally, their Pliicker matrices are given by 

L = M(p3,P2-|-(Pl)*, L = M(p3,P2-(Pl)* 

where M and the * operator were defined in section ISTI 
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